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THE CHARGING PROCESS 
IN A HIGH- SPEED , S INGLE- CYLINDER , FOUR- STROKE ENGINE 
By Blake Reynolds, Harry Schecter, and E. S. Taylor 

SUMMARY 



Experimental measurements and theoretical calcula- 
tions nave been made on an aircraft- type , s i Ug 1 e- cy 1 i n& er 
engine, in order to determine the physical nature of the 
inlet process, especially at high piston speeds. The en- 
gine was run at speeds from 1,500 to 2,600 r.p.m. (mean 
piston speeds 1,370 to 2,380 feet per minute). Measure- 
ments were made of the cylinder pressure during the inlet 
stroke, and of the power output and volumetric efficiency, 
xvieasur e.nent s were also made, with the engine not running, 
to determine the resistance and mass of the air in the in- 
let-valve port at various crank angles. 

A theoretical analysis of the process was made, on a 
purely mechanical basis (no account being taken of thermal 
effects). The results of this analysis indicate that, as 
far as the mechanics is concerned, mass has an appreciable 
effect but friction plays the major part in restricting 
flow . 

Prom the pressure records it was found that trie cylin 
ier was filled to very nearly atmospheric pressure at bot- 
tom center, even at the highest speed. (Similar results 
were noted on a larger cylinder running at piston speeds 
up to 2,750 feet per minute.) The observed fact that the 
volumetric efficiency is considerably less than 100 per- 
cent is therefore attributed to thermal effects. An esti- 
mate was made of the magnitude of these effects in the 
present case, and their general nature discussed. 

In connection with the analysis, it was necessary to 
obtain accurate records of the absolute cylinder pressure. 
Modifications of the ri.I.T. indicator for this purpose are 
described in an appendix. 
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INTRODUCTION 



The maximum power output of an engine is closely re- 
lated to the quantity of air which it can pump. At high 
piston speeds, this quantity may he limited "by the flow 
characteristics of the inlet valve. Previous experimenters, 
recognizing this limitation, have devoted considerable ef- 
fort to determine the characteristics of steady flow (ref- 
erences 1, 2) and dynamic flow (references 3, 4, 5, 6) 
through valves. Work has also "been done (references 3, 7, 
8, 9) on the effect of inlet pipes of various lengths on 
the dynamics of air flow. This work has "been largely de- 
voted to engines operating at lower piston speeds, and in- 
formation relating to the flow in the valve itself under 
actual operating conditions is incomplete. 

The purpose of the present investigation is to estab- 
lish the factors governing the flow of air* into the cyl- 
inder of an internal combustion engine operating at high 
piston speed, with particular attention to the relative 
importance of mass and friction in controlling the flow. 
It was expected that results of such a study might very 
well lead to a better understanding of the inlet process, 
if not to modifications of the conventional inlet arrange- 
ment for the purpose of increasing air capacity. 

APPARATUS AND EXPERIMENTAL PROCEDURE 



The engine used was an experimental Pratt and Whitney 
Wasp, j r , air-cooled aircraft engine cylinder mounted on a 
universal single-cylinder crankcase (fig. 1). (Dimensions 
are given in appendix I.) Cooling air was drawn (not 
blown) over the cylinder, to avoid the possibility of pre- 
heated air entering the inlet port. 

Indicator diagrams were taken by means of an improved 
!!«.!. T« balanced-pressure indicator, using a Famboro type 
cylinder unit developed especially for these tests. (See 
appendix II.) 

During the tests the inlet pipe was removed, the fuel 



♦For convenience, the word "air" is used to denote the 
fuel-air mixture in the inlet system. 
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being injected directly into the inlet port by means of a. 
nozzle located as shown in figure 1 and supplied from a 
special high-speed Bosch fuel- in j ect ion pump. Spark tim- 
ing and' fuel-air ratio were adjusted for best power con- 
ditions, The exhaust pipe wa.s also removed during the. 
tests, since it was found that the exhaust system caused 
considerable hack pressure on the engine (see fig. 2). 

Figure 3 shows' a set of records of cylinder pressure, 
taken at various speeds from 1,500 to 2,600 r.p.m. These 
pressure curves are traced from light-spring indicator 
diagrams which were taken with a pointer speed of approxi- 
mately 0.1 inch per second. The indicator contained S.A.3. 
No. 20 oil, and the corresponding error due to indicator 
piston friction is about 0.03 inch or 0.15 pound per square 
inch. (See appendix II, and figure 23.) 

Quantities required for the theoretical analysis were 
determined as follows. 

The natural frequency, u) n /2TT, of the cylinder and 
valve acting as a Helmholtz resonator was determined at 
various points over the inlet stroke by exciting the reso- 
nator by" means of a loudspeaker placed near the valve port 
and actuated by a calibrated beat-frequency oscillator 
(General Radio type 713-3). The results are accurate to 
within a few cycles per second, as resonance was quite 
sharp and easily discernible. Figure 4 shows the variation 
of the natural frequency with crank angle. 

The effective valve area A defined as the area of 
an equivalent rounded-edge orifice, was determined from, 
steady state flow data. The cylinder was evacuated to 
maintain a pressure drop across the intake valve and port 
of 12 inches of water by an exhauster connected to the ex- 
haust port, the exhaust valve being held open. The quan- 
tity of air passed per unit time was measured by an air 
meter in the exhaust line. The meter was an K.A.C.A. 
Roots supercharger with a capacity of 0.150 cubic foot per 
revolution. Figure 5 shows the variation of the lift and 
of the effective area with crank angle. 



Volurnetr i c Ef f i ci ency 



Volumetric efficiency e y is herein defined as the 
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ratio of the mass of air taken in during one suction stroke, 
to the .-ciass of air which would fill the piston displacement 
at inlet pressure and temperature.* For convenience in the 
discussion, the term "volumetric deficiency, 11 defined as 
(1 - e v ), will also be used. 

Because of the difficulty of making a direct measure- 
ment of volumetric efficiency without disturbing trie phe- 
noaiena under investigation, it was decided to compute the 
volumetric Efficiency from power measurements. 

With the inlet pipe connected, the engine was operated 
at the same speeds as in the tests with no inlet pipe. A 
Hash "HyTor 11 compressor, which delivered air at approxi- 
mately atmospheric temperature, was used to increase inlet 
pressure until the power was approximately the same as had 
been recorded with the inlet pipe off. Prom air measure- 
ments made by" means of an 1; . A . 'C . A . Hoots supercharger used 
as a meter, it was found that the volumetric efficiency 
was very nearly equal to 0.0053 (i.m.e.p.), where the 
i. m.e.p. was determined by adding the brake m.e.p. to the 
motoring friction m.e.p. The volumetric efficiencies used 
in this report are all computed from the above relation, 
and should be accurate to about 2 percent. 



T3330BETICAX ANALYSIS 

According to the simplest view, the velocity of the 
air in the inlet valve is proportional to the piston ve- 
locity, the proportionality coefficient "being the ratio of 
piston area to effective valve area. Such a view has been 
used in a previous investigation (reference 9), the results 
of which show that this gives a good approximation to the 
actual motion. A acre careful examination requires that 
trie impedance to flow at the valve be taken into account, 
flora the flow may be impeded by the inertia of tho air in 
the valve, and by the friction encountered by the air in 
going through the valve. 

For the purpose of analyzing the motion, the air in 
the valve is considered to move as a solid plug, of mass 
ill. This plu^r is acted upon by three forces: 



*For the experimental conditions used, inlet pressure and 
temperature are atmospheric. 



[ 
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B* , the force due to the pressure outside the cyl- 
inder , 

F o , the force due to the pressure inside the cyl- 
inder , 

!F 3 , the force due to the resistance to flow through 
the valve. 

If u denotes the velocity of flow through the valve 
(measured positive for flow into the cylinder), the equa- 
tion of motion is 



the forces being measured positive in the direction of u. 

The effective mass m may he found, approximately, 
from engine dimensions and the measurements described on 
page 3. The stiffness of the air in the cylinder is (ref- 
erence 10, pp. 201-2C2). 

k = ^ (2) 

where 7 is ratio of specific heats Cp / 0 V 

F, pressure in cylinder 

■ A, effective valve area 

V, volume of cylinder 

The angular frequency of the free vibrations of the air in 
the valve is 



iD n =V*/ai (3) 



Thus m may be obtained from tiie relation ffi = k/ju n 2 . The 

value of k is calculated from equation (2), using atmos- 
pheric 7 and p since the frequency measurements were 
made at atmospheric conditions; and the value of (, -> n c? is 

obtained from figure 4. The resulting curve showing the 
variation of m with 6 is given in figure 5. 
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The force due to the pressure outside the cylinder 
is simply 



(4) 



where F 0 is atmospheric pressure. 

The force due to the pressure inside the cylinder is 



ff a = -AP 



(5) 



In order to determine r it is assumed that the expansion 
of the air in trie cylinder takes place i sentropi cally . 
The relation between pressure and volume, for the cylinder 
contents at any instant, is then 



4? _ V *Y 
F V 



Audt 



(6) 



since the change in volume of the cylinder contents is 
(dV-Audt) and the volume is V. Integrating, and solving 
for P , 



P 9 P 



J c (7Ac/v)dt 



(7) 



where P Q is initial pressure, assumed atmospheric, 

V n , initial volume, assumed equal to clearance 
volume 

The force due to the resistance to flow through the 
valve may be approximated "by 



F, - -A(AP) 



(a) 



where AP is the pressure difference (measured positive 
in direction of u) required to overcome tho resistance 
and maintain steady flow of velocity u. For low veloci- 
ties of flow this pressure drop is equal to ^PU a (where 
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P is the density of the air in the valve) , Under actual 
running conditions, however, the velocity reaches quite 
lar a*e values, SO that it is desirable to use the more ex- 
act expression. This expression may readily be derived if 
it is assumed that the expansion of the air from atmospher- 
ic conditions Up to the conditions in the valve takes place 
i sen t r opi cal ly and in accordance with the perfect-^as la^s. 
The relation is (reference 11, p, 225) 

a? = p 0 {i ~ (i - i- ?v i §a u 2 ) 7TT } (9 ) 

where P Q is atmospheric density. 

fjT convenience, the equation of motion i 1 1 "be written 
in terms of a velocity r a tip (instead of the velocity), 
namely 

$ = u/c (10) 

wher e 

- V?o/N (ID 



Making appropriate substitutions, the equation of motion 
is obtained in the form 



it .n dt P + mc v 1 J " vTV 

(12) 



or 



d 3 

ae 



/ (Q/6J0pae 

If + 1 - Y - Z f r ' 



(13) 



where 6 is crank an^le ( :l c gr ees A . T , C . ) 
N, engine speed (r.p.m.) 



ra dc 
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X = AP 0 /rac 

7-1 

Y = 1 - fl - 

\ 2 

Z x (V /V) 7 
~~ 3 



p 8 ) 



Q = YAc/v 

The equation representing the motion is a nonlinear 
differential-integral equation with variable coefficients. 
Since the coefficients cannot readily be expressed analyt- 
ically, it is practically impossible to solve the equation 
except by graphical-numerical methods. This can be done 
rather simply, since the equation gives the slope (dp /& 9) 
in terms of p and 6. 

In order to solve the equation of motion (13), it is 
necessary to determine how the coefficients M, X, Z, 
and Q vary with 0 # from the solution of the equation, 
the pressure in the cylinder may be obtained simply by us- 
ing equation (7) . 



Determination of Coefficients 

The coefficient M may readily be determined from m 
(e.g., graphically, by plotting log m against 0 and 
measuring the slope of the resulting curve). The variation 
of si with 9 is shown in figure 6, 

Tne coefficient X may readily be obtained from the 
data in figure 5, The variation of X with 8 is shown 
in figure 7. 

The function Y may readily be calculated. A bino- 
mial expansion is convenient, three terms being sufficient 
(i.e., up to p 6 ) « Tnis function is plotted in figure 8, 
for three different values of 7. 

The coefficient Z is readily obtained from engine 
dimensions; this coefficient is plotted in figure 9. 

The coefficient Q,, plotted in figure 10, is readily 
obtained from engine dimensions and the data in figure 5. 

Careful examination reveals that the Mj3 term belongs 
in the equation of motion only when M is positive. When 
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M is positive the plug of air is picking up mass from 
rest , and the force which accelerates this masi in a part 
of the force, acting on the plug, so that the |£g term must 
be retained in tne equation. On the other hand, when 1.1 
is negative , the mass lost by the plug is acted upon by a 
force, out tuis force is not a part of the force on the 
plug, and is irrelevant to the motion of the plug. There- 
fore when U is negative, th.e Mj3 terra does not "belong 
in tho equation of motion. 

A simplified analogy is the case of a moving train 
picking up or discharging bars of sand. If the moving 
train has bags thrown onto # it (from rest), the train will 
bo decelerated (a > 0, urn ~ -mu) ; but if the moving 
train has^bags thrown off it, the train will not be accel- 
erated (m < 0, mx = 0) . 

The whole argument, of course, hinges on tne supposi- 
tion that the moving mass discarded by the plug (when 
id < 0) is brought to r^st by external forces -which act on 
the discarded mass without affecting the remaining mass 
(i.e., the plug). This will not always be true; e.g., if 
there i s a. sizable inlet pipe, the moving mass will be 
brought to rest partly by the containing vails adjacent to 
the valve, and partly by mixing with the air along the 
walls of the pipe; this latter part will affect the pros- 
sure acting on tne plug. In such cases the effect is *orop~ 
erly expressed by retaining some portion of the US tern, 

In the experimental work herein described, however, 
the inlet oipe was completely removed. Since the length 
of the port is only a few inches, tne L!j3 term was omitted 
when A < 0. 

figure 11 (lotted line) is a pressure curve resulting 
from the solution of equation (13) for 2,500 r.p.m., rith 
the experimental pressure record superimposed for compari- 
son. It appears that the pressure is much too low except 
at tne beginning of the suction stroke. This led to an in- 
vestigation of nonme chani cal effects on volumetric effl- 
c i ency . 

THEHliAL EFFECTS 

Exarninat ion of the experimental indicator diagrams re- 
vealed the rather astonishing fact that tne press ure at the 
time of inlet c 1 o s i n ™ s n g a r 1 y c on s t an t over tho r a n g e of 
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speeds investigated. In order to check this observation, 
the pressure Pg was measured at a number of points- along 

the early part of the compression stroke. From each of 
these pressures, an equivalent pressure at bottom center 

P was obtained by means of the relation PqVq 1 * 35 = 

p V 1,35 where V, is the cylinder volume at bottom cen- 

11 1 
ter and Vg is the cylinder volume at the crank angle 9 

for which Pg was measured. The values obtained by this 
method were consistent with each other, indicating that 
the exponent 1.35 was properly chosen for the cylinder 
contents at this time. The average of the equivalent F ± 
is plotted in figure 12. (It should be noted that the ac- 
tual pressure at bottom center is not P 1 since the inlet 

valve is still ojjen and t lie flow lg continuing after bot- 
tom center.) Results of this measurement show that p, 
is nearly 14.7 pounds per square inch and changes little 
over the speed range investigated. Further confirmation 
of this fact was found by noting tnat the compression lines 
of diagrams taken at various speeds appeared to be identi- 
cal when superimposed, A series of three indicator dia- 
grams taiien from a Wright Cyclone engine with inlet pipe 
removed, at 1,600, 2,000, and 2,400 r.p.m. (piston speeds 
1,830, 2,290, 2,750 feet per minute), showed a similar 
lack of variation in P 1 . 

The question immediately arises as to what causes the 
volumetric deficiency and what makes it vary with r.p.m. 
It is clear that the cause must be either an increase in 
temperature as the air comes into the cylinder, or a re- 
duction of the partial pressure of fresh air in the cylin- 
der due to dilution with fuel or residual exhaust gas. 

Effect of Residual Gas 

Part of the volumetric deficiency could be due to the 
effect of (1) the residual gas in the cleara nee space Do- 
ing at a pressure higher than atmospheric (3) the fact that 
the specific heat of the residual ,::as is somewhat different 
from that of the fresh charge, which will result in a dif- 
ferent volume before and after mixing the two gases. These 
effects have been investigated (appendix III) and found to 
be negligible for the cases under consideration. 
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Effect of Fuel Evaporation 

The addition of fuel also has an effect on volumetric 
efficiency. If no fuel evaporation occurs the effect is 
negligible. Evaporation causes a decrease in partial pres- 
sure of air and a decrease in the temperature of tho mix- 
ture. The latter effect is much larger. Complete evapor- 
ation of fuel under adiabatic conditions and at the fuel- 
air ratio used (0.080) would cause a net volumetric defi- 
ciency of -0.04. (See reference 12.) 



Effect of Throttling on the Temperature T x 

Assuming, for the moment, adiabatic conditions: as 
air flows through the inlet valve, the temperature will be 
reduced "by expansion. The velocity it acouires during this 
process will., however, "be completely converted into turbu- 
lent flow and eventually into random molecular motion. As- 
suming that the conversion of energy of turbulent flow to 
heat energy is completed immediately after entry into the 
cylinder, the temperature of the entering air will be the 
same as the outside temperature since air may be considered 
as a perfect gas for the pressures and temperatures under 
consideration (Joul e- Thomson experiment). The air which 
has entered will then follow the pressure changes which oc- 
cur in the cylinder and will eventually arrive at a temper- 
ature corresponding to isentropic compression from the en- 
tering pressure to the cylinder pressure at, say, bottom 
center. Thus, in general, each element of air will arrive 
at a different temperature, and since the pressure at bot- 
tom center is very nearly atmospheric, the temperature of 
most of the fresh charge in the cylinder will be above at- 
mospheric. An estimate of the average temperature ri so 
due to this process is possible if we take the piston work 
during the suction stroke equal to the heat energy causing 
the average rise in temperature. 

The measured mean effective pressure during the suc- 
tion stroke at 2,500 r.p.m. is 2.5 pounds per square inch 
below atmospheric and the work' done is therefore 0.465 
B.t.u. per cubic foot of inlet air. With an inlet-air den- 
sity of 0.052 pound per cubic foot, 5 percent residual 
gas, and a mean specific heat of 0.24 B.t.u. per pound, 

the average temperature rise is P* = 30° F. 

• '0.0S2 X 1.05 X 0.24 

This will account for a volumetric deficiency of about 

0.05. It is questionable how much of this effect is pres- 



12 



Itf.A.C. A. Tec;ndcal Kote JTo.. 575 



ent in the actual case since a considerable proportion of 
the energy of flow may not have been converted into random 
energy when the inlet valve closes. 

Effect of Heat Transfer from the Cylinder Walls 

There remains but one possible explanation for the 
observed volumetric deficiency: There must be a consider- 
able transfer of heat to the charge, during t}xe suction : 
stroke, from the surfaces of the cylinder, inlet port, and 
valve. No data are available showing how much heat is trans 
ferred during flow through the port and. valve and how much 
heat flows to the charge in the cylinder, nor is any exper- 
imental information. Available as to how the rate of heat 
flow varies during the cycle. It seems' reasonable to sup- 
pose that the rate of heat flow from the inlet valve and 
valve seat to the flowing gas is comparatively lar^e be- 
cause of the- high relative velocity, and large temperature 
difference. At the beginning of the suction stroke the 
heat flow within the cylinder is undoubtedly in a direction 
from the gas to the cylinder wall, but as fresh charge 
comes in the te nrp 6 r a ture of the cylinder contents will be 
rapidly reduced and neat will flow in the opposite direc- 
t i on . 

It should be noted that the effect of heat transfer 
may vary with speed and therefore be partly responsible 
for the variation in volumetric efficiency with speed. In 
addition to the fact that the nature of flow is Changing, 
the temperature of the sv.rfaces.of the cylinder walls and 
particularly the temperature of the inlet valve may vary 
considerably with speed. 

CORRECTION 0? ANALYSIS TO • ALLOW ?0H THERMAL EFFECTS 



•An attempt was made to take account of the various 
thermal factors influencing flow. In order to accomplish 
this*, it was necessary to find how much volumetric defi- 
ciency was due to effects other than a change in total 
pressure at the time of inlet closing. This was done by 
multiplying the observed volumetric efficiency by the ob- 
served ratio 14.7/P 1# This quantity has been plotted in 
figure 12. At 2,500 r.p.ra. the pressure is responsible for 
only a small volumetric deficiency and 14.7 0- fl? 1 = 0.78, 

Assuming that the flow of incoming air is adiabatic through 
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the valve and that immediately thereafter the net effect 
of heat transfer, fuel evaporation, and kinetic energy- 
transfer is to raise its absolute temperature, thereby re- 
ducing its density oy a constant f ac t or , " 0 . 78 , the ecuation 
of motion is changed only in that the integral in equation 
(13) is divided by 0,78. The solid line in figure 11 is 
the pressure computed using this correction. The agreement 
with the experimental record is much improved, indicating 
that the above assumptions are reasonable. ■ 

Results of similar computations made for 2,200 and 
1,500 r.p.m. are plotted, with the exper i mental records su- 
perimposed for comparison, in figures (13) and (14). The 
equivalent pressure P 2 obtained from the computed pres- 
sure at 198° crank angle, is plotted in figure 12. The 
difference between the computed and measured P- is ap- 
proximately constant' and equal to 0.3 pound per square inch. 
Apparently the analysis checks the variation of ? 1 with 
speed over the range investigated. In this connection, it 
is interesting to note that the initial" solution using no 
correction for thermal effects is approximately equivalent 
to an increase of 28 percent in r f p.m. and indicates a con- 
siderable reduction in P at this higher speed (approxi- 
mately 3,300 r.p.m.). 

previous investigators (ref ereneas Ji 5) ^ av e assumed or 
concluded that the friction plays the major part in restrict- 
ing the flow through the valve, the mass being relatively 
unimportant. In order to determine just how much the iner- 
tia does or does not affect the inlet process, the effect 
of omitting the mass from the foregoing analysis has been 
tried. With zero mass, the equation of motion becomes 

- • / (Q/6H)pd9 

1 - Y - Z e (14) 



The resulting curve of pressure is plotted in figure 15 
with the curve from figure 11 and the experimental record. 
It is evident that the analysis including mass fits the ex- 
perimental data better, although the two solutions seem to 
be nearly identical during the latter part of the stroke. 
The nature of the difference between the curves suggests 
that an intermediate analysis, omitting the term dp/d9 but 
retaining the term HQ , might be satisfactory. 
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CONCLUSIONS 

The following conclusions have "been drawn from the 
measurements and analysis: 

1. pressure drop in the valve is responsible for an 
extremely small part of the volumetric deficiency on the 
cylinder tested, even at 2,380 feet per minute mean piston 
speed . 

2. The rise in temperature of the charge coming 
through the inlet valve and in the cylinder before inlet 
closing is responsible for a relatively large volumetric 
deficiency, 

3. It is possible to calculate the pressure drop 
in the valve with good accuracy, 

4. The friction of the air in the valve port is re- 
sponsible for most of the pressure drop. 



Re cornraendat ions 

It is recommended that further study on this problem 
be directed particularly toward a study of the heat trans- 
fer process. 

It also seems advisable to continue the work to still 
higher piston speeds. (This was not possible in the pres- 
ent investigation on account of the valve gear, which fi- 
nally failed at 2,600 r.p.m.) The actual valve-lift curve 
should be observed since the valve lift when the engine is 
running at high speed may be considerably different from 
the lift measured statically at the same crank angle (ref- 
erence 13) • 

Further improvement in the method of measuring pres- 
sures is desirable, 

Massachusetts Institute of Technology, 
Cambridge, Mass., June 1938. 
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APPSSPIX I 

NOTATION AND NUMERICAL DATA 
Notation 

A, effective inlet valve area 

e, velocity of sound, •AP 0 /P 0 [ 

q , "volumetric efficiency 

k, stiffness of air in cylinder 

L, connect in^-roi length 

rn, effective mass of air in the inlet valve 

M - I 45 

ni d9 

N, engine speed, r ,p # m. 

P, pressure in cylinder during suction stroke 

P , atmospheric "ores sure 
o 

P , equivalent pressure in cylinder at bo t ton) center 

F , pressure in cylinder at D during compression stroke 
9 

Q « 7Ac/v 

r, crank radius 

R , compression ratio 

t , t i rae • 

T , t emperatur e 

T. , temperature of cylinder contents at bottom center 

u, velocity of air in the inlet valve 

V , volume of cylinder during suction stroke 

V , volume of cylinder at bottom center 
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clearance volume 
displacement volume 

volume of cylinder at 8 during compression stroke 
AP Q /mc 



7 




(V 2 /V) r 
u/c 

base of natural logarithms 
ratio of specific heats, 0 /C y 
crank angle, degrees A.T.C. 
density of air in valve 
atmospheric density 

natural frequency of air in valve and cylinder, t/kjm 
physical Constant s 

= 530 °H. 

= 14.7 l'b . / sq . in . 

The average value of 7 for the cylinder contents 
(air and residual gas) during the suction stroke 
has been taken as 1,3. For the charge at the end 
of trie suction stroke 1.35 has been used. 

1 , 130 f t . /sec . 
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Engine Dimensions 

Sore 5-3/16 inches 

Stroke 5-1/2 inche s 

Crank radius r 

■ ~ 0.2'34 

Connecting rod length 

Piston displacement, V d ...... 113.7 CU, in. 

Compression ratio, R . , , 5.5 

Inlet-valve seat, inner diameter .... 2 inches 

From the geometry, 

- 1 _ - cos e + S in 2 e \ 
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APPENDIX II 
M.I.T. BALANCED -PRE S SURE INDICATOR 



The indicator used to obtain pressure versus crank- 
angle diagrams of the cylinder was an improved Ivl.I.T. type 
"balanced-pressure indicator (reference 14), manufactured 
by the Otico Instrument Company (figs. 16, 21). The baU 
anced-pressure principle upon which thi s - indi cat or works 
is shown diagramma t i cal ly in figure 17. Typical light- 
and heavy-spring diagrams are shown in figures 18-20. 

Since the determination of the small pressure drop 
across the inlet valve required an accurate record of the . 
absolute pressure in the cylinder, greater precision was 
called for than had previously been accomplished. There- 
fore, considerable experimental work was done on the in- 
dicator in an attempt to increase its accuracy and deter- 
mine its range of error. 

As far as the indicator itself was concerned, it was 
found that somewhat greater accuracy could be realized 
for light-spring diagrams if the pressure was not in 

operation during the time that the diagram was being re- 
corded, since the oressure fluctuations from the pump were 
enough to broaden the recorded line. Therefore, the pump 
was used only to drive the pointer to its maximum pressure 
position and was then shut off, the card being taken as the 
spark pointer returned to its atmospheric position. In 
taking heavy-spring diagrams, the width of the line is not 
increased perceptibly by the action of the pump, so that 
it is not necessary to use this lengthier procedure. 

A second source of inaccuracy in the indicator itself 
arises from the speed at which the spark pointer is allowed 
to return to its atmospheric position. The indicator pis- 
ton, spring, and oil-filled cylinder together constitute 
an overdamped vibratory system which has a definite time 
lag in following any forced motion. The constants of this 
system were determined by a simple experiment, as follows: 
The indicator drum was turned at a convenient speed, and 
the piston was moved away from its atmospheric position by 
hand, the bleed valve being open and the pump off. The 
"atmospheric" switch was then pressed and the piston re- 
leased. The "atmospheric" switch of the Otico model (fig. 
21) produces 120 sparks per second from the spark point to 
the drum so that the position of the pointer was recorded 
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every 120th .of a second, The rotation of- the drum spread 
the spark holes so that they could "be easily identified 
and their distance from the starting position measured. 
A curve of displacement versus time (fig. 22) was then 
plotted, from which the time constant of the system was 
determined (equation (3), below) . 

As' a forcing function a constant rate of change of 
"balancing pressure was assumed. This assumption approxi- 
mates fairly closely the actual operating conditions, 
where the "bleed valve is opened continually wider as the 
balancing pressure decreases. From the assumed forcing 
function, the error was calculated as a function of the 
pointer speed and the spring constant. The error is ex- 
pressed in inches displacement, on the diagram, from the 
position corresponding to the actual "balancing pressure. 
Figure 23 is a plot of the error against the reciprocal of 
the speed of the spark pointer, for various springs with 
oil in the cylinder and for the 5-pound spring with and 
without oil. As may he seen, the removal of the oil raises 
the permissible speed considerably, in the case of light- 
spring recording. 

The type of pick— up unit which had previously given 
the best results is one utilizing the motion of a clamped 
diaphragm (fig. 24) to make and break the grid circuit of 
the thyratron. The diaphragm, stamped from spring steel 
stock two- to f iv e —thousandths inch thick, is supported on 
either side by perforated steel disk?, so that it will with- 
stand a pressure somewhat more than 1,000 pounds per square 
inch. It was found by experimenting with several different 
designs that the nonad jus table contact-disk feature with a 
flat surface (fig. 24) makos for maximum ease of cleaning 
and reassembly without the need of readjustment. All of 
the required clearance for diaphragm motion is supplied in 
the front disk, which is cut down two -thousand ths of an inch 
from a supporting rim at the edge. The back disk is flat. 
Approximately sixty holes, no. 55 drill, with no grooving, 
makes for good construction and operation without intro- 
ducing any discernible error. The possibility of error 
arises from the fact that when the diaphragm is pressed 
against one of the disks the effective area may not be the 
same on both sides of the diaphragm. The possible effect 
of this unbalance, however, was not noticeable in a care- 
ful check between units with grooved and ungrooved disks. 
Therefore, it has been assumed negligible. ■ 

The chief disadvantage of the d iaphragm- type pick-up 
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unit is that it is subject to a zero error, i.e., the en- 
tire indicator diagram is shifted up or down with respect 
to the atmospheric line. This error arises from the fact 
that the diaphragm . general ly has an initial pressure 
against the contact disk or against the front disk; a cer- 
tain pressure is required on one side or the other in or- 
der to place the diaphragm in such a position that a slight 
change in pressure will make or "break contact. The pres- 
sure required to place the diaphragm in this position is 
the zero error of the unit and is equal to the displace- 
ment of the diagram with respect to the atmospheric line. 

In an attempt to measure the magnitude of the zero 
error so that corrections for it could be made, several 
means of calibration wore devised. The simplest, an appa- 
ratus for static calibration (fig. 26), gives only a rough 
measure of the error, say to within a half inch of mercury. 
The shortcoming of the method lies apparently in the fact 
that the make and break occur much more slowly than under 
running conditions, a gradual change of resistance causing 
an indeterminate time of spark. As a means of increasing 
the rapidity of make and break, the fluctuations of the in- 
dicator pump can be utilized as follows. The unit is at- 
tached to the indicator in the usual manner except that the 
cylinder side of the diaphragm is exposed to atmospheric 
pressure. ffith the pump off and the bleed valve open, the 
drum is rotated and the atmospheric line recorded; this 
line will be in the correct position, corresponding to 
equal pressures on both sides of the unit. Then the bleed 
valve is closed and the balancing pressure slowly varied, 
until the unit interrupts the current and records an at- 
mospheric line; this line will be in an incorrect position, 
corresponding to the pressures on the two sides of the 
unit differing by the zero error. Thus the difference in 
position between the two linos is the zero error of the 
unit. Because of the pressure fluctuations, the width of 
the atmospheric line recorded by a well-adjusted unit may 
be as much as 0.5 pound per square inch, but an average 
line drawn through this band was found to be in error by 
not more than 0.25 pound per square inch as compared with 
the reversal method of calibration. 

Since neither of the foregoing methods simulates run- 
ning conditions very closely, it was doomed wise to intro- 
duce another method, namely the reversal method, which 
would do so. One cylinder of an Indian motorcycle engine 
was redesigned with a two-to-one compression ratio head so 
that a complete indicator diagram could be obtained on a 
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scale sufficiently large to insnre good accuracy. llo 
valves were used, the pressure level being controlled "by 
the flow of air into the cylinder through a s:nall capil- 
lary tube. A critically damped manometer, although it 
did not read the average pressure in the cylinder correct- 
ly, nevertheless did serve as an excellent means of indi- 
cating -whether the average pr e s ?ur e r emai nod constant. 

. Tne method consisted of taking two diagrams on one 
card: one with tne unit connected in the conventional man- 
ner and the other with the unit reversed end for end, i.e., 
the hack side of the diaphragm, which is usually exposed 
to "balancing pressure, was connected to. the cylinder, while 
tne front side was connected to the pressure line. Thus, 
whereas in one case tne zero ^rror wculd raise the diagram 
with respect to the atmospheric line, in the other case it 
would lower the diagram by the same amount, Tne resulting 
difference between the two diagrams is twice the zero error 

Since the flow of air in the unit is restricted by a 
small tube on the indicator side of the diaphragm, whereas 
there is no restriction cn the cy linior s ide (fig. 24), it 
was necessary to design an adapter for the latter side, so 
as to make the- assembly symmetrical about the diaphragm. 
The adaptor, of courrc, introduced considerable lag in the 
transmission of pressure variations, so that true cylinder 
conditions were no longer present at the diaphragm. The 
indicator, however, records whatever variations are present 
at the diaphragm, and since those variations are unchanged 
by -the reversal of a symme t r i cal assemoly, the "forward" 
and "reverse 11 diagrams from the unit are congruent. The 
calibration is unaffected by the lag. 

In practice it was found that the two diagrams were 
not exactly congruent due to the fact that the adapter did 
not make the unit quite symmetrical. Therefore, the zero 
error was taken to be half tne difference between the aver- 
age pressure lines of the two diagrams, A typical reversal 
calibration is shown in figure 27. 

As these methods of calibration seemed to give fairly 
consistent results among themselves , i was decided to in- 
vestigate tne effect of temperature on tne zero error. The 
results of the investigation proved conclusively that the 
diaphragm units under considerate on we-e unsuitable for the 
measurement of absolute pressure to a greater accuracy than 
two pounds per square inch. The thought of accurate cali- 
bration of this fcype (i.e., the . clampedr-diaphragm type) was 
therefore , abandoned . 
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As an alternative, the freely floating valve (Farn- 
boro type) unit was resorted to (fig. 25). In order to 
obtain a clear-cut card from this type of unit it was 
found necessary to attach a fine wire between the valve 
and the center terminal; the contact between the valve 
stem and the necessarily loose guide was apparently poor 
enough to cause a broadening of the lines of the diagram. 
Even with this wire connection, however, the reversal 
method showed no zero error. It is reasonable to assume 
that a charge of temperature will not introduce a zero 
error in this type of unit. This assumption is certainly 
justified by the consistency of the results obtained from 
this type of unit. 

Another possible source of error is introduced with 
the valve unit, namely, a time lag due to the inertia of 
the valve. However, no perceptible lag was found, either 
from a careful check between diagrams taken with these 
and the diaphragm units, or from a study of the symmetry 
about the top-center line of a diagram from a nonfiring 
engine. The entire diagram was, of course, taken on 
"break" so that the time for the valve to cross the gap 
would not affect the recording. 

The valve unit is, then,, an excellent method of deter- 
mining absolute pressures and as such can be used to deter- 
mine the variation in the zero error of a diaphragm unit 
under operating conditions. The procedure is to superim- 
pose diagrams taken by the two typos of units under the 
same operating conditions. The true atmospheric line, de- 
termined from the diagram given by the valve unit, can 
then be drawn on the diagram from the diaphragm unit and 
the zero error measured off directly. By using this method 
over a range of operating conditions in the Wasp engine, it 
was found that the zero error of a diaphragm unit may change 
as much as 1 # C pound per square inch with a variation in 
speed from 1,500 to 2,500 r.p.m.. A test previously run on 
a C.3P.B. engine showed an increase in diaphragm temperature 
from 425° to 500° F, as the speed was increased from 900 to 
1,800 r.p.m.; therefore, it is possible that the observed 
change in zero error is duo to the effect of temperature. 

A corrroar i son of the results of the various methods of 
calibration is shown in table I. From those and other 
data it is safe to draw the following conclusions on cali- 
brat ion. 

1) The static calibration - either hot or cold - of 
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a diaphragm unit that has been "run in u can be 
relied upon conservatively for an accuracy 
within 5 pounds per square inch - in most cases 
- within 3 pound? per square inch. 

2) Calibration with the indicator pump is comparable 

to static calibration, 

3) The reversal method of calibration does not war- 

rant the additional time required (i.e., not 
for routine calibration). 

4) The valve unit needs no calibration, and can 

therefore be used to calibrate the diaphragm 
unit as mentioned above. 

The effect of vibrations on the width of the recorded 
lines was studied in two ways. The C.F.R. engine with its 
heavy water-jacketed cylinder was considered substantially 
free of vibration, particularly when running at a low com- 
pression ratio. Diagrams taken from it, then, were as- 
sumed to be the best that could be obtained from the par- 
ticular units which were used. In few cases did either 
the valve or diaphragm units give lines on the light-spring 
diagram which were wider than 0.3 pound per square inch, 
except during the exhaust stroke where it is doubtful 
whether the cycle repeats closely. The same units in the 
ffasp cylinder under the worst conditions of vibration gave 
lines sometimes as broad as 1.0 pound per square inch, but 
usually on the order of 0.7 pound per square inch. The 
valve units showed a tendency toward slightly broader 
lines than did the diaphragm type. Where lines broader 
than 1,0 pound per square incii were encountered, it was 
found that cleaning the unit would reduce the width to 0,7 
pound per square inch, or loss, for normal running condi- 
tions of the Wasp cylinder. 

The second method of testing for the effect of vibra- 
tion consisted of recording the make or break of contact 
of the unit when placed in a blind plug mounted in the 
cylinder wall. This shielded the unit from cylinder pres- 
sure but subjected it to vibration from the cylinder, 
ffith the engine running, the balancing pressure was slowly 
reduced to atmospheric. When the pressure was close to 
atmospheric in the case of the valve unit or close to the 
zero-error pressure in the case of the diaphragm unit, the 
motion of the vibrating cylinder head forced the unit away 
from the valve or diaphragm, the latter acting as seismic 
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masses. Thus contact was made or "broken and a "band was 
recorded on the diagram. These hands were compared with 
similar hands taken with the units stationary. Again the 
broadening effect of the vibration was found to be on the 
order of 0.5 pound per square inch. 

The effect of vibration on the units, then, is merely 
to broaden the indicator lines. An average line drawn 
through the broader one is probably as accurate as the 
finer line from a nonvibrat ing engine. 

TABLE I 

INDICATOR UNIT CALIBRATION 

All pressures - Inches of Mercury 

Static Reversal 

Jo Id Hot Cold 



Unit Diaphragm 
No. thickness 
(in.) 



Running 
(against 
valve 
uni t ) 



3 


0.005 


13 


9 


10 


10.5 


7.3 


5 


• .002 


11 


11 


12 


12.3 


9 


7 


.002 


. 5 


.4 


.5 


1.0 


.4 


2 


.005 


+ 1 . 5 


- .5 


+ 1.5 


+ 1.6 


-4.0 


6 


.002 


5 


3 


2 


2.5 


3.5 




Un i t 


No. 5 


, diaphragm 0. 


005 inch 






Calibrated 


agai 


nst valve 


uni t 


in Wasp 


engine 



r . p . m . 

1 , 500 

1 , 800 
2, 000 

2 , 200 
2,400 
2, 600 



Inches Hg , 



7 

25 
0 
6 
5 

25 



Unit No 



Static 



o 
6 

10 
10 

3 
6 

5 



1.5 
1.7 

2.3 

5.0 

20 



Pump 



Rev er sal 



1.3 
1.4 
5.8 

■ .8 

2.8 
20.3 



6.9 
2.8 
.5 
5.0 

20.4 
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Calculation of Error Due to Pointer Speed 

The error due to the pointer speed may "be calculated 
"by considering the response of the pointer to a forcing 
function. 

The equation of motion is 

mx 4- cx 4- kx =s kXf (1) 

where x is pointer position 

m, effective mass of pointer, spring, etc. 

0 | damping constant 

k, spring constant 

X> , forcing function 

Since the system is overdamped, equation (1) may "be 
written in the approximate form 

i six f (2) 

T T 

wher e 

T= c/k 

The time constant, t, may be determined from figure 22, 
which shows the measured response of the system upon be- 
ing released from an initial ..deflection x Q . The comple- 
mentary equation representing the free motion of the sys- 
t em i s 



x + i x = 0 



solving for 1/t, 



i * I dx _ d In x 

T £ = ~ M dt dt 



Thus 1/t is simply proportional to the slope of the 
straight portion of the response curve, figure 22. 
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: ,7ith the time constant known-, it is possible to find 
the error due to trie lag in the response of the moving 
system to the varying "balance pressure which acts as the 
driving force. As a close approximation to operating con- 
ditions, it is assumed that the pointer moves at a constan 
speed v. The forcing function is thus 

x f - vt (4) 

With this forcing function, the steady-state solution* of 
the equation of motion is 

x * v(t - t) (5) 

The correct position of the pointer, corresponding 
to the actual value of the balancing pressure, is X* 

(i.e., the position which the pointer Would have if the 
system were purely stiffness-controlled; or, with the ac- 
tual system, the position it would have if the motion 
were infinitely slow). The recorded position is x, giv- 
en by equation (5). The error € is defined as the cor- 
rect position minus the recorded position, 



From equations (4) and (5), 

* =• v T (6) 



The spring constant k is found** from the spring 
rating I (pounds per square inch per inch), 

k m SK (7) 

where S is the area of the indicator piston. 

* Since the time constant is of the order of 0.1 second, 
the complementary solution becomes negligible within about 
0.5 second after motion starts. 

**0r it can be measured directly, if necessary. 
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The damping constant c ( = kf) is presumably inde- 
pendent of the spring constant, and therefore may he de- 
termined by using any spring, knowing k and measuring 
T as described above. Once this has been done, the valuo 
of t for any other spring may readily be found from tho 
relation t = c /k by using the appropriate value of k. 

From measurements made using the 5-pound spring, with 
S.A.E. No. 20 oil in the indicator cylinder, the error was 
found to be 

» * 1.37 I 

K 

where e is in inches, v in inches per second, and K 
in poinds per square inch per Inch, A similar measurement, 
without oil, gave a result of <€ = 0.11 v for the 5-pound 
spring. 

These results are plotted in figure 23. 



APPENDIX III 

CALCULATION OF THE EFFECT OF RESIDUAL EXHAUST GAS 
ON VOLUMETHIC EFFICIENCY 



An estimation of the effect of residual exhaust gas 
on volumetric efficiency may be made for adiabatic condi- 
tions, that is, no heat flow to the walls of the container, 
and under the assumption that the residual exhaust gases 
are i son tropi cally compressed or expanded to the same 
pressure as the inlet gases and then mixed. 

The following symbols will be used: 

P 0 , atmospheric pressure 

T 0 , atmospheric temperature 

P , pressure in cylinder at end of suction stroke (bottom 
center ) 

T , temperature of mixture at end of suction stroke 
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P , pressure in cylinder at ond of exhaust stroke (top 
center ) 

T , te^erature of residual gas at end of exhaust stroke 

6 ' 

7 , cylinder volume at end of suction stroke 



V o , clearance volume 

R - f,/fg i compression ratio 



n r 



number of mo Is of air. in cylinder at end of suction 
stroke 

number of mols of residual gas in cylinder at end of 
suction stroke 



G , mean molal heat capacity of air (constant pressure) 



■a 



C , mean molal heat capacity of residual gases (constant 



r 



pressure ) 



'7, the ratio of specific heat at constant pressure to 
specific heat at constant volume, for residual 
gas at temperature T 6 

R Q , universal gas constant 

When the volumetric efficiency is unity, n %9 the number 
of mols of air taken in will be 



V V 1 - V 2> 



*0 T o 



The volumetric efficiency e is therefore 



(1) 



From the gas law, 

? V = n r R n T. (2) 

6 2 r- 0 6 

P V = (n a + n r ) R 0 I (3) 



N.A.C.. A. Technical LTote Ho. 575 



31 



If the residual gag is expanded or compressed isen- 
tropically to the pressure Pj and the temperature T 7 
and then mixed with the incoming charge, 

1 - 7 



% - t 6 ( u ) y (4) 



The temperature of the. mixture I may be found from the 
r elat ion 



or 



C r n r 



From (2) and (3) , 



- = E 1 (6) 

p 6 T i 



Combining (4), (5), and (6) , 



T o G r * P. J Me 
~ = 1 2 (7) 

1 1 a r 6 1 

R - 



CoaMning (o) with (1) and (2) 

? T 



e v = -i 

P 0 B - 1 



and substituting (7) in (8), 
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e v - — (9) 

"o R - 1 



It will "be noted that one of the unknowns T 1 re- 
mains in the expression. However, T g is not altogether 
independent of I. f and T 1 /f fi may "be considered con- 
stant without appreciable loss of precision since T 1 /T 6 
is small, (0 /C - 1) is small, and the product is small 

compared to the other terms. The equation in this form 
is considerably simpler than it would he if T x were re- 
placed by its value in terms of T Q • 

Sample Computat ions 

The ratio of the mean molal heat capacity of residual 
gas to that of air is approximately 1.1. For the tempera- 
ture range investigated T % /T 6 is approximately 0.3 and 
7 is about 1.31, Substituting these values in (9), 

PR- 1.1(? /? ) ***** + G.0Z(P 6 /P 1 ) 

e = -1 S 1 - (10) 

P 0 R - 1 

or, for the engine used (R ~ 6.5), 
p 

e v - — [1.18 - 0.20(P 6 /P 1 )°' 76 * 0.005 ^/^) ] (11) 
? o 

In order to estimate the minimum effect of residual 
gas, the following values wore determined from an indicator 
diagram, taken at 2,000 r.p.m., which showed a minimum val- 
ue of P . 

6 

-~ = 1 0.965 
o r o r i 



Substituting these values gives 
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e v = 1.13 - 0.20 (0.965) 0 - 76 + 0.005 



"V 

= 0.99 



A second set of values was determined from an indi- 
cator diagram, taken at 1,500 r.p.m., which showed a max- 



imum value of P 6 . 



P, ? P 

— =1 -1 = JL = 1.027 

0 0 1 



The result of substitution of tnese values is 

e v = 0.985 



It was concluded that the effect of residual gas 
could be neglected. 
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Fig. 1 




FIGURE l.-WASP JR. CYLINDER ON UNIVERSAL CRANKCASE 
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Figure 2.- Indicator diagrams showing effect of exhaust pipe on cylinder pressure. 8200 r.p.m. 
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Figure 3.- Records of cylinder pressure. 1200 to 2F00 r.p.m. 
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Figure d.- Variation of valve lift, effective valve area, and effective mass, with crank angle. 
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Figure 11.- Comparison of indicator diagram with calculated pressure, showing effect of 
correcting analysis for thermal effects. 3600 r.p.nw 
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Figure 12.- Variation of Pi and e y witn r.p.m. 




Figure 13.- Comparison of indicator diagram with calculated pressure. (Corrected for 
14.7e v /Pi= 0.80.) 2200 r.p.m. 
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Figure 14,- Comparison of indicator diagram with calculated pressure. (Corrected for 
I4.7e y /P x * 0.81.) 1500 r.p.m. 




Figure 15,- Comparison of indicator 
mass. (Corrected for 14 



diagram 



with calculated pressure, showing effect of neglecting 
- 0.78.) 2600 r.p.m. 
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Fig. 16 




FIGURE 16. -M. I. T. BALANCED PRESSURE INDICATOR 




Figure 17.- Schematic drawing of indicator. 




Figure 18.- 

Wasp, jr. 1800 R.P.M. 
(With exhaust pipe) 
Diaphragm type indicator unit 



Pressure scale 5 lb «/ 3< l* ln « 

inch 



Figure 19.- 

C.P.R. 1200 H.P.M. 
Valve type indicator unit 



Pressure scale 5 lb «/ 3 <l» 

inch 
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ELEMENTARY wiring diagram - 
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Figure 21 
Wiring diagram 
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WIRING DIAGRAM 
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produce slight pressure 
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Figure 26.- Static calibrator. 



CO 



Reversal Calibration Method 



Forward 



Ave Press. For. 




Atmos. R 



\04 



^Reverse 



Area Under Forward » 9.45 in 2 

Average Pressure * *£n*h * above base 

line 

Area Under Revier^e - 9 . 40 in* 
Averacje Pressure * \,04" 

E= 2 * Zero error ^ 1.02" 

Figure 27. 



Press. Scale- 5 lb/in* = \ in. 
Zero Error ^. * 5.2 in. H9. 



